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Abstract 

A multidimensional extremal problem in the idempotent algebra 
setting is considered which consists in minimizing a nonlinear func- 
tional defined on a finite-dimensional semimodule over an idempotent 
semifield. The problem integrates two other known problems by com- 
bining their objective functions into one general function and includes 
these problems as particular cases. A new solution approach is pro- 
posed based on the analysis of linear inequalities and spectral proper- 
ties of matrices. The approach offers a comprehensive solution to the 
problem in a closed form that involves performing simple matrix and 
vector operations in terms of idempotent algebra and provides a ba- 
sis for the development of efficient computational algorithms and their 
software implementation. 
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1 Introduction 

The development of solution methods and computational algorithms for 
solving multidimensional extremal problems is one of the important con- 
cerns in the linear tropical (idempotent) algebra [H El O HI El [7] . The 
problems under consideration involve the minimization of linear and non- 
linear functionals denned on finite-dimensional semimodules over idempo- 
tent semifields and may have additional constraints imposed on the feasi- 
ble solution set in the form of linear tropical equalities and inequalities. 
Among these problems are idempotent analogues of the linear programming 
problems [SI El [7] and their extensions with nonlinear objective functions 
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[ini HD H21 US HH US H6l HZl HE]. There are solutions to certain prob- 
lems where both objective function and constraints appear to be nonlinear 

da [2D]. 

Many extremal problems are formulated and solved only in terms of one 
idempotent semifield, say the classical semifield K m ax,+ in 0[2[T7]. Some 
other problems including those considered in [HJ E21 EE EH E3 El ES] are 
treated in a more general setting, which includes the semifield R m ax,+ a s a 
particular case. Furthermore, proposed solutions frequently take (see, e.g. 
[TU1 El [201 El EZ] ) the form of an iterative algorithm that produces a solution 
if any, or indicates that there is no solutions, otherwise. In other cases, as 
in [HI [101 [T2J [131 EH E3 [13 IS], direct solutions are given in closed form. 
Finally, note that most of the existing approaches offer some particular 
solutions rather than give comprehensive solutions to the problems. 

In this paper, we consider a multidimensional extremal problem that 
is a generalization of the problems examined in [10l HI E3 [HI 13 (13 ES] ■ 
Particular cases of the problem arise in various applications including growth 
rate estimation for the state vector in stochastic dynamic systems with event 
synchronization \10\ [P3] and single facility location with Chebyshev and 
rectilinear metrics [15} 116] . On the basis of implementation and further 
development of methods and techniques proposed in [TQ1 ED E21 E3 EH E3 
[TBI [TBI I2T] , we give a complete solution to the problem in a closed form that 
provides an appropriate basis for both formal analysis and development of 
efficient computational procedures. 

The rest of the paper is as follows. We begin with an introduction to 
idempotent algebra and outline basic results that underlie the subsequent 
solutions. Furthermore, examples of tropical extremal problems are pre- 
sented and their solutions are briefly discussed. A new extremal problem 
is then introduced and a closed-form solution to the problem under general 
conditions is established. Finally, solutions to some particular cases and 
extensions of the problem are given. 

2 Basic Definitions and Results 

We begin with basic algebraic definitions and preliminary results from 
[T2"l [T3] to provide a formal framework for subsequent analysis and solutions 
presented in the paper. Additional details and thorough investigation of the 
theory can be found in [JJ El El H El El [7] . 

2.1 Idempotent Semifield 

We consider a set X that is closed under binary operations, addition © and 
multiplication (g), and equipped with their related neutral elements, zero 
and identity 1. We suppose that the algebraic system (X, 0,1,©, ®) is a 
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commutative semiring with idempotent addition and invertible multiplica- 
tion. Since for all x £ X + , where X + = X \ {0}, there exists its multiplica- 
tive inverse x~ l , the semiring is commonly referred to as the idempotent 
semifield. 

The integer power is introduced in the usual way to represent iterated 
multiplication. Moreover, we assume that the power with rational exponent 
is also defined and so consider the semifield to be radicable. 

In what follows we omit the multiplication sign ® as it is usual in the 
conventional algebra. The power notation is always used in the above sense. 

The idempotent addition naturally induces a partial order on the semi- 
field. Furthermore, we assume that the partial order can be completed to 
a total order, thus allowing the semifield to be linearly ordered. In the fol- 
lowing, the relation signs and the min symbol are thought of as in terms of 
this linear order. 

Examples of linearly ordered radicable idempotent semifields include 

Kmax,+ = (K U {-oo}, -oo,0, max, +), 
Kmin,+ = (K U {+00}, +oo,0, min, +), 
K m ax,x = (K+U{0},0,l,max, x), 
Kmin,x = (K+ U {+00}, +00, l,min, x), 
where R is the set of reals, R + = {x € R\x > 0}. 

2.2 Idempotent Semimodule 

Consider the Cartesian product X ra with column vectors as its elements. A 
vector with all components equal to is called the zero vector and denoted 
by . The operations of vector addition © and scalar multiplication <g> are 
routinely defined component-wise through the scalar operations introduced 
on X. 

The set X n with these operations forms a finite-dimensional idempotent 
semimodule over X . 

A vector is called regular if it has no zero components. The set of all 
regular vectors of order n over X + is denoted by X™ . 

For any nonzero column vector a; = (ij) £ X" we define a row vector 
x = (x~) , where x~ = x^ 1 if X{ 7^ , and x7 = <D otherwise, i = 1, . . . , n. 

2.3 Matrix Algebra 

For conforming matrices with entries from X, addition and multiplication 
of matrices together with multiplication by scalars follow the conventional 
rules using the scalar operations defined on X . 

A matrix with all entries that are equal to is called the zero matrix 
and denoted by . A matrix is row (column) regular if it has no zero rows 
(columns). 
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Consider the set of square matrices X nxn . As in the conventional alge- 
bra, a matrix is diagonal if its off-diagonal entries are equal to . A diagonal 
matrix that has only 1 as the diagonal entries is the identity matrix and 
denoted by /. Finally, the exponent notation stands for repeated multipli- 
cation for any square matrix A with the obvious condition that A = I. 

The set X nxn with the matrix addition and multiplication forms an 
idempotent semiring with identity. 

For any matrix A = (a^) its trace is given by 



tr.4 = 0a, 



i=l 



A matrix is called reducible if simultaneous permutations of rows and 
columns put it into a block-triangular normal form, and irreducible other- 
wise. The normal form of a matrix A € X nxn is given by 



A 



( A X1 



\ 



(1) 



\ A sl A s2 ... A ss ) 



where An is either irreducible or zero matrix of order n 
is an arbitrary matrix of size rij x nj for all i = 1, . . . , 
n\ -\ \-n s = n. 



whereas Aij 



j < i. 



and 



2.4 Spectrum of Matrices 

Any matrix A 6 X nxra defines on the semimodule X n a linear operator with 
certain spectral properties. Specifically, if the matrix A is irreducible, it has 
a unique eigenvalue that is given by 

n 

A= Qtr 1 /™^" 1 ) (2) 

m=l 

whereas all corresponding eigenvectors are regular. 

Let the matrix A be reducible and have the form (JTJ) . All eigenvalues 
of A are among the eigenvalues Aj of the diagonal blocks An , i = 1, . . . , s . 
The value A = Ai © • • • X s is always an eigenvalue, it is calculated by ([2]) 
and called the spectral radius of A . 



2.5 Linear Inequalities 

Suppose there are a matrix A £ x mxn and a regular vector d £ X™ . The 
problem is to solve with respect to the unknown vector x £ X™ the linear 
inequality 

Ax < d. (3) 
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Clearly, if A = (D, then any vector a; is a solution. Assume now the 
matrix A ^ to have zero columns. It is easy to see that each zero column in 
A allows the corresponding element of the solution vector x to take arbitrary 
values. The other elements can be found from a reduced inequality with a 
matrix that is formed by omitting zero columns from A and so becomes 
column-regular. The solution to the inequality for column-regular matrices 
is as follows. 

Lemma 1. A vector x is a solution of inequality ([3|) with a column-regular 
matrix A and a regular vector d if and only if 

x < (d~A)~. 

For a given square matrix A S x nxn and a vector 6 G X", we now find 
all regular solutions x £ to the inequality 



Ax © b < x. 



(4) 



To solve the problem, we follow an approach based on the implemen- 
tation of the function Tr that maps each square matrix A 6 X nxn to a 
scalar 



For any matrix A £ 



Tr(A) = 0tr^ m . 

m=l 

" x "' we introduce a matrix 



A* = I © A © ■ ■ ■ © A 



n-l 



Assume a matrix A to be represented in its normal form ([TJ. We define 
a diagonal matrix 

/ A u 

D = 



and a low triangular matrix 



A 



21 



A 



s,s— 1 



\ 



which present the diagonal and triangular parts of the decomposition of A 
in the form 

A = D®T. (5) 



Note that if the matrix A is irreducible, we put D = A and T 
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Theorem 1. Let x be the general regular solution of inequality (JU) with a 
matrix A in the form of ([5]) . Then the following statements are valid: 



1. If Ti{A) < 1, then x = (D*T)*D*u for all ueX" such that u>b. 

2. If Tr(.A) > 1, then there is no regular solution. 

3 Tropical Extremal Problems 

We now turn to the discussion of multidimensional extremal problems for- 
mulated in terms of idempotent algebra. The problems are established to 
minimize both linear and nonlinear functionals defined on semimodules over 
idempotent semifields, subject to constraints in the form of linear equalities 
and inequalities. 

In this section, the symbols A and C stand for given matrices, b, d, 
p, q, g and h for vectors, and r and s for numbers. We start with an 
idempotent analogue of linear programming problems examined in [9, 7J 
and defined in terms of the semifield K m ax,+ to find the solution x to the 
problem 



A solution technique that is based on an iterative algorithm and called 
the alternating method is proposed which produces a solution if any, or 
indicates that there is no solution otherwise. 

The technique is extended in [IT] to provide an iterative computational 
scheme for a problem with nonlinear objective function given by 



There are certain problems which can be solved directly in a closed form. 
Specifically, an explicit formula is proposed in [8] within the framework of 
optimization of max-separable functions under disjunctive constraints for 
the solution of the problem 



Furthermore, in \10\ fT2j [131 IT^l I15| . a problem is examined which is to 
find regular solutions x that provide 




Ax®b <Cx®d. 



min (p x(Br)(q x © s) 
Ax b < Cx © d. 



-i 



T 

mm p x, 
Cx > b, 
g < x < h. 



min x Ax. 
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To get a closed form solution to the problem, an approach is applied that 
uses results of the spectral theory of linear operators in idempotent algebra. 

Finally, a closed form solution based on a technique of solving linear 
equations and inequalities is derived in [18] for the problem 

min (x~pQ)q~x), 
Ax < x. 

In the rest of the paper, we consider a problem with a general objec- 
tive function that actually contains the objective functions of the last two 
problems as particular cases. For the problem when there are no additional 
constraints imposed on the solution, a general solution is given in a closed 
form. 

4 A New General Extremal Problem 

Given a matrix A 6 x nxn and vectors p, q G X n , consider the problem to 
find x that provides 

min (x~ Ax x~p © q x). (6) 

A complete explicit solution to the problem under general conditions as 
well as to some particular cases and extensions is given in the subsequent 
sections. 

4.1 The Main Result 

We start with a solution to the problem in a general setting that is appro- 
priate for many applications. 

Theorem 2. Suppose A is a matrix in the form ([T]), p is a vector, q is a 
regular vector, A is the spectral radius of A, and 

A=(q~p) 1/2 , / u = A©A/0. 

Define a matrix 

A^ = fi A = T^, 

where and T /t are respective diagonal and lower triangular parts of A^ , 
and a matrix 

b = (d^td;. 

Then the minimum in (0) is equal to [i and attained if and only if 

x = Bu 

for all regular vectors u such that 

fj, p < u < fi(q~B)~. 
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Proof. We show that both A and A are lower bounds for the objective 
function in ©, and then get all regular vectors x that yield the value 
fi = A © A of the function. To verify that A is a lower bound, we write 

x ~ Ax © x~p © q x > x~ Ax. 

Assume the matrix A to be irreducible and A to be its unique eigenvalue. 
We take a corresponding eigenvector xq and note that for all x € X™ , 

XqXq = 1, XXq > (x~ Xq)^ 1 I . 

Furthermore, we have 

X Ax = X~ AxXq Xq = X~ A(xXq~ )xq 

> x~Axq(x ^Xq)" 1 = \x~xq(x ~Xq) _1 = A. 

Consider an arbitrary matrix A taking the form (pQ). Any vector x 
now admits a decomposition into subvectors xi,...,x a according to the 
decomposition of A into column blocks. With the above result for irreducible 
matrices, we obtain 

si s s 

x Ax = CD x ■ A-ijXj > CD x ■ A{{X{ > CD ^ = 

i=l j=l i=l i=l 

Now we show that A = (q^p) 1 ^ 2 is also a lower bound for the objective 
function. We have 

x~ Ax © x~p © q~x > x p © q x. 

Let us take any vector x £ X" and denote 

r = x~p © q x. 

The last equality leads to two inequalities 

r > q x > 0, r > x~p. 

Multiplication of the first inequality by r~ 1 x~ from the right gives x > 
r q xx > r q~ ■ Substitution of x~ > r~ l q~ into the second results 
in r > r ' q p = r _1 A 2 , whence it follows that 

x~p © q~ x = r > A. 

By combining both bounds, we conclude that 

x Ax © x~p © q x > A © A = ^. 
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It remains to find all regular solutions x of the equation 

x Ax © x p © q x = {A. 

Since x~ Ax(Bx~ p@q x > fi for all x G X™ , the set of regular solutions 
of the equation coincides with that of the inequality 

x Ax © x p © q x < [i, 

which itself is equivalent to the system of inequalities 

x~Ax © x~p < n, (7) 
q x < \i. (8) 

Let us consider inequality ([7|). After multiplication of the inequality by 
fi~ 1 x from the left, we write 

A^x © [i~ x p < pT x xx~ Ax © yT x xx~p < x, 

and then arrive at the inequality 

A^x © fi^p < x. 

On the other hand, left multiplication of the obtained inequality by jix~ 
directly yields inequality ([7]), and thus both inequalities are equivalent. 

Since Tr^) = Tr(/z -1 A) < 1, we can apply Theorem [1] to the last 
inequality so as to get the general solution of inequality ([7]) in the form 

x = (d;t^d; u = Bu, 

where u € X" is any vector such that u > [i~ 1 p. 

Substitution of the solution into inequality flHJ) gives an inequality q~ Bu < 
[i. Application of Lemma [1] to the last inequality yields u < n(q~ B)~ . 

By combining lower and upper bounds obtained for the vector u, we 
finally arrive at the solution 

x = Bu, 

for all u € X™ such that 

H _1 p<u< n(q~B)~, □ 

4.2 Particular Cases and Extensions 

Consider problem Q with an irreducible matrix A. Since in this case 
Dp = Afj, , = (D , and B = A* , the statement of Theorem [2] takes a 
reduced form. 
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Corollary 1. If A is an irreducible matrix, then the solution set of © is 
given by 

x — -A.. u> 

for all u € X™ such that 

^p^uK^q-AlY. 

Specifically, when A = 0, we have B = A*^ = I and fj, = A, The 
solution set is further reduced to 

A l p < x < Aq, 

which coincides with that in |18| . 

Suppose the vector q in problem (|6|) is irregular. In this case, the matrix 
q B in the inequality 

q~Bu < /j, 

may be not column-regular, which prevents direct application of Lemma [1] 
as in Theorem [2J 

Let J = supp(q~.B) be the set of indices of nonzero elements in the 
row vector q~B. Denote by (q~B)j and uj subvectors that have only 
components with indices from J. The solution of the above inequality is 
given by the constraints uj < fi(q~B)j for the subvector uj , whereas the 
rest components of the vector u can take arbitrary values. 

Now we can somewhat weaken conditions of Theorem [2] as follows. 

Theorem 3. Under the assumptions of Theorem^ let q ^ be arbitrary 
vector and J = supp(q~ B) . 

Then the minimum in Q is equal to fi and attained if and only if 

x = Bu 

for all regular vectors u such that 

lT l p<u, uj<[i(q~B)j. 

Finally note, that when q = we have J = and so the upper bound 
for u disappears. 

5 Conclusion 

A complete closed-form solution has been derived for a tropical extremal 
problem with nonlinear objective function and without constraints. The 
solution actually involves performing simple matrix and vector operations 
in terms of idempotent algebra and provides a basis for the development of 
efficient computational algorithms and their software implementation. 
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As a suggested line of further research, solutions to the problems under 
constraints in the form of tropical linear equalities and inequalities are to 
be considered. Practical examples of successful application of the results 
obtained are also of great interest. 
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